Measurement of the vector and tensor analyzing powers for dp- elastic scattering at 

880 MeV 



P.K.Kurilkin a , V.P.Ladygin a , T.Uesaka b , K.Suda c , Yu.V.Gurchin 3 , A.Yu.Isupov a , K.Itoh d , M.Janek ae , J.-T.Karachuk af , 
T.Kawabata b , A.N.Khrenov a , A.S.Kiselev 3 , V.A.Kizka 3 , V.A.Krasnov ag , N.B.Ladygina 3 , A.N.Livanov a ' g , Y.Maeda h , 
A.I.Malakhov a , S.M.Piyadin 3 , S.G.Reznikov 3 , S.Sakaguchi b , H.Sakai', Y.Sasamoto b , K.Sekiguchi c , M.A.Shikhalev a , 

T.A.Vasiliev a , H.WitalaJ 



o 

(N 



in 



"VBLHEP-JINR, 141980 Dubna, Moscow region, Russia 
b Center for Nuclear Study, University of Tokyo, Tokyo 113-0033, Japan 
c RIKEN Nishina Center, Wako, Saitama 351-0198, Japan 
''Department of Physics, Saitama University, Saitama, Japan 
e Physics Department, University of Zilina, 010 26 Zilina, Slovakia 
f Advanced Research Institute for Electrical Engineering, Bucharest, Romania 
g Institute for Nuclear Research, Moscow, Russia 
h Kyushi University, 6-10-1 Hakozaki, Higashi-ku, Fukuoka-shi 812, Japan 
'University of Tokyo, Bunkyo, Tokyo 113-0033, Japan 
'Institute of Physics, Jagiellonian University, Krakow, Poland 



Abstract 



1 The vector A y and tensor analyzing powers A yv and A xx for dp- elastic scattering were measured at T . = 880 MeV over the cm. 



rlab 

qj[ angular range from 60° to 140° at the JINR Nuclotron. The data are compared with predictions of different theoretical models based 
.on the use of nucleon-nucleon forces only. The observed discrepancies of the measured analyzing powers from the calculations 
CJ require the consideration of additional mechanisms. 
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1. Introduction 

The nucleon-deuteron (Nd) elastic scattering is the simplest 
composite-particle scattering process. It has been widely used 
to test how well the few-nucleon systems can be described in 
terms of nucleon-nucleon (NN) interactions. At energies be- 
low the n production threshold (E=210 MeV/u), Faddeev cal- 
culations provide rigorous descriptions of the scattering pro- 
cess. During the last several years, polarization observables 
in Nd elastic scattering have been studied in a number of ex- 
periments at RIKEN CI H S Hi KVI Jl S S S S E d, 
IUCF HI Hdl and RCNP Til [Till. The considerable ex- 
perimental activities were stimulated, in particular, by the ob- 
served discrepancy of ~30% between differential cross section 
data |Hl[l3l measured at energies of 65-135 MeV/nucleon and 
results of Faddeev calculations fljll . Many of the discrepancies 
for the differential cross sections and vector analyzing pow- 
ers at these energies are remedied by the inclusion of the 2n 
-exchange three-nucleon forces (3NFs) such as TM-3NF fl^l. 
UrbanaIX-3NF S or TM99 jSj] into the calculations H. 
But theoretical calculations with 3NFs still have difficulties in 
reproducing data of some spin observables, for instance, tensor 
analyzing powers. 

At high energies, the forward angle Nd elastic scatte ring 
has been successfully described by the Glauber approach 12311 



which takes into account both single and double scattering 
terms. The interference between the single- and double- scat- 
tering amplitudes including the D- state in the deuteron wave 
function (DWF) allowed one to explain the filling of the cross 
section diffractive minimum 1 2411 . 



At energies 200-600 MeV/nucleon ; however, large discrep- 
ancies between the experimental data 11151 1 16ll25l 12611 and theo- 
retical predictions in the minimum of the differential cross sec- 
tion are persistent even after inclusion of 3NFs [22j,|27|]. There 
are several possible origins for these discrepancies. One is rel- 
ativistic effects: the relativistic Faddeev approach |28, 29| has 
been developed for intermediate energies range. It was found 
that when only NN forces are included in the calculations, the 
relativistic effects are significant only at backward scattering 
angles. They are relatively small in the minimum of the differ- 
ential cross section where the discrepancies between 2N force 
predictions and data are largest B28I1 . The extension of the Fad- 
deev calculation technique into the relativistic regime [ 3(il 31 ] 
also does not provide a reasonable description of the experi- 
mental data, only with NN forces. Recently, however, it is re- 
ported that large changes of the elastic scattering cross section 
in the region of angles ranging from the minimum of the cross 
section up to very backward angles are observed when 3NFs 
are included in the calculations 1291. 
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Another possibility is the manifestation of new 3NFs. Exist- 
ing 3NF models 1 19, 2(1 21] mostly deal with low-momentum 
diagrams, while diagrams with higher momentum can make 
visible contributions in scatterings at 200-600 MeV/nucleon 
where large momentum transfers are relevant. It is also pos- 
sible to consider reaction mechanisms which are not included 
in the present Faddeev calculations. In the energy region con- 
sidered, effects due to n production and on-shell A excitation 
are the candidates. 

Thus it is strongly anticipated that a reaction theory for the 
nucleon-deuteron scattering at 200-600 MeV is established and 
further information on new 3NF is extracted from the scattering 
data. The energy region lies above the n- production threshold 
and the applicability of the Faddeev calculation techniques, at 
least as it is, is not trivial. On the one hand, these energies are 
not large enough to apply the Glauber approach. 

In this work new results on the analyzing powers A y , A vy and 
A xx for dp- elastic scattering measured at Tf°= 880 MeV over 
the cm. angular range from 60° to 140° are presented and are 
compared with several theoretical calculations, including Fad- 
deev [18] and Glauber M23I1 ones. 



2. Experimental procedure 

The experiment has been performed at the Internal Target 
Station (ITS) 13211 at the superconducting synchrotron, Nu- 
clotron at the Laboratory of High Energy Physics, Joint Insti- 
tute for Nuclear Research. The ITS consists of a spherical scat- 
tering chamber and a target sweeping system. The scattering 
chamber is fixed on the flanges of the Nuclotron ion tube. The 
disk mounting six different targets is located on the axle of a 
stepper motor. A target used for the measurement is moved to 
the center of the ion tube when the particles are accelerated up 
to the required energy. A 10 pm CH2 film was used as a pro- 
ton target. A carbon wire was used to evaluate the background 
originating from the carbon content in CH2. The signal from 
the target position monitor f33ll was used to tune the accelera- 
tor parameters to bring the interaction point close to the center 
of the ITS chamber. The signals from the monitor were also 
stored as raw data so that one can use the position information 
in off-line analysis. 

The polarized deuteron beam was provided by the atomic- 
beam polarized ion source POLARIS [34]. Nuclear polariza- 
tion is provided via radio-frequency (RF) hyperfine transitions. 
In this experiment the data were taken for three spin modes: un- 
polarized, "2-6" and "3-5", which have theoretical maximum 
polarizations of (pz,Pzz) = (0,0), (1/3,1) and (1/3,-1), re- 
spectively. Two different RF cells of POLARIS [34] with the 
working frequencies of 384.9 MHz and 320.1 MHz have been 
used to provide the "2-6" and "3-5" transitions, respectively. 
The spin modes were cycled by spill-by-spill. The polarized 
deuteron beam was accelerated up to T l ^ b -%^,Q MeV by keep- 
ing the quantization axis perpendicular to the beam-circulation 
plane of the Nuclotron. The typical beam intensity in the Nu- 
clotron ring was 2-3 xlO 7 deuterons per spill with a duration of 
~ 1 s, independently on the spin mode. 



The detection system was designed for analyzing powers 
measurements in a wide range of initial deuteron energies [35]. 
The detector support with mounted 46 plastic scintillation coun- 
ters was placed downstream the ITS spherical chamber. Each 
plastic scintillation counter was coupled to a photo-multiplier 
tube Hamamatsu H7416MOD. Nine proton detectors were in- 
stalled for left, right and up, but due to space limitation - only 
four for down directions. The proton detectors were placed 
at a distance of 600 mm from the target. The angular span 
of one proton detector was 2° in the laboratory frame, which 
corresponds to ~4° in the cm. Four deuteron detectors were 
placed at scattering angles of deuterons coinciding kinemati- 
cally with the protons for left, right and up scattering. Only one 
deuteron detector was used to cover the solid angle correspond- 
ing to down scattering. In addition, one pair of detectors was 
placed to register two protons from quasi-elastic pp- scattering 
at 9 PP =90° in the cm. in the horizontal plane to monitor the 
beam luminosity. The deuteron- and quasi-elastic pp- detectors 
were placed at a distance of 560 mm from the target in front of 
the proton detectors. 

The scattered deuterons and recoil protons were detected in 
kinematical coincidence over the center-of-mass angular range 
of 60-140°. The analyzing powers A y , A yy and A xx were mea- 
sured at nine (eight) different angles in the cm. defined by the 
position of the proton counters placed in the horizontal (verti- 
cal) plane. 
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Figure 1 : The time difference between the signals for deuteron and proton de- 
tectors at 7j*=880 MeV and a scattering angle of 70° in the cm. a) - the 
normalized spectra obtained for CH2 and carbon targets (open and shaded his- 
tograms, respectively); b) - the result of the CH2-C subtraction. The dashed 
lines are the prompt timing windows to select the dp- elastic scattering events. 



Selection of the dp- elastic scattering events was based on 
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the kinematical coincidence of the scattered deuterons and re- 
coil protons by the scintillation counters at several angles in the 
cm. The energy loss correlation and time-of-fiight difference 
for the signals from the corresponding proton-deuteron detector 
pairs were used. The time difference between the signals from 
the deuteron and proton detectors for CH2 and carbon targets 
gated by the correlation on the energy losses are shown in the 
upper panel in Fig.[T]by the open and shadowed histograms, re- 
spectively. The dashed lines in Fig.[T]represent the prompt tim- 
ing window for the selection of the dp- elastic scattering events. 
The carbon contribution is mostly due to quasi-free scattering 
and appears as a broad bump. It is eliminated partly by the 
criteria on the energy losses correlation in deuteron and proton 
scintillation detectors. However, the contribution from the car- 
bon content in CH2 inside the timing window varied between 
10 and 25% depending on the detection angle. The final selec- 
tion of the dp- elastic events has been performed by the CH2-C 
subtraction of the normalized time-of-fiight difference spectra 
for each pair of the conjugated deuteron and proton detectors 
for each spin mode of the PIS H34I1 . The quality of the CH2-C 
subtraction is demonstrated in the bottom panel in Fig.[T] 

The beam polarization has been measured before the data 
taking using the asymmetry of the dp- elastic scattering yields 
at Tj b =270 MeV with the same detection system [36]. Al- 
though the CH2 film was used as a proton target, no measure- 
ments with a carbon target were made, because the background 
from the carbon content in CH2 was found negligibly small 
1 35], being in good agreement with RIKEN measurements [Q]]. 
Data at several scattering angles were used to decrease the sta- 
tistical errors of the beam polarization. The values of the an- 
alyzing powers Ay, A yv , A xx and A xz at these angles were ob- 
tained by the cubic spline interpolation of the precise RIKEN 
data Ell]. The values of the vector p y and tensor p yy com- 
ponents of the beam polarization for the "2-6" and "3-5" spin 
modes of POLARIS H34fl are presented in Table Q] The system- 
atical errors indicated in Table Q] are due to the uncertainty of 
the dp - elastic scattering analyzing powers at T l ? b =270 MeV 



^Pzz{A yy {8)cos 2 (t> + A xx (0)sin 2 <f>)], 
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Table 1 : The values of the vector p y and tensor p vv beam polarizations for the 
"2-6" and "3-5" spin modes of POLARIS [ 34] obtained at T l f=270 MeV (H. 



Spin 
mode 


Py 


Apf 


A P r 


Pyy 


Apyf 




"2-6" 
"3-5" 


0.216 
0.208 


0.014 
0.011 


0.002 
0.002 


0.605 
-0.575 


0.024 
0.020 


0.005 
0.005 



The analyzing powers A y , A yy and A xx for the dp- elastic scat- 
tering at 7^ fl6 =880 MeV were measured simultaneously. The 
analyzing powers are defined in terms of the xyz coordinate of 
z\\ki, y\\ki xkf, and x\\y x z, where k t and kf are the incident and 
scattered deuteron momenta, respectively. The yields of the lip- 
elastic scattering with the vector pz and tensor pzz deuteron po- 
larizations can be written as 13811 



where N po i(6, (p) and Nq(0, <p) are the yields corrected by the 
beam luminosity and dead-time with polarized and unpolarized 
beams, respectively, 6 is the scattering angle, and <p is the az- 
imuthal angle with respect to the beam direction. The azimuthal 
angles <p for the detectors placed in the directions of left, right, 
up, and down are 0, n, -n/2, and n/2 radians, respectively. 
The analyzing powers A y , A yy and A xx were extracted by us- 
ing the normalized yields n{8,<p) = N po i{8, (f>)/No{8, (/>), defined 
in Eq. © 



A y {8) 
A yy (8) 
A xx (0) 



n(0, 0) - n(0, n) 
3pz 

n(6, 0) + n(0, n) - 2 
Pzz 

n(8, -n/2) + n(6,n/2)-2 
Pzz 



(2) 



N po i(e,4>) 



No(6,4>)-[l + -pzA y (e)cos(t 



Such a method does not require an accurate knowledge of the 
detector geometries and/or efficiencies of the detection system. 
The analyzing powers A y , A vv , and A xx extracted by Eq. (0 were 
averaged over "2-5" and "3-6" spin modes of POLARIS Ir34ll . 

3. Results and discussion 

The angular dependencies of the vector A y and tensor ana- 
lyzing powers A vv and A xx obtained at T' d ab =8d,0 MeV are pre- 
sented in Fig. [2] by the solid symbols. The error bars are the 
statistical only. The systematic errors for the vector A y and ten- 
sor analyzing powers A yv and A xx from the uncertainty of the 
normalization of the beam polarization are ~ 7% and ~ 4%, 
respectively. The systematic errors due to the procedure of the 
dp- elastic events selection were estimated to be ~ 5.5% and 
~ 2.5% for vector and tensor analyzing powers, respectively. 

The solid lines in Fig. [2] are the results of the nonrelativistic 
Faddeev calculations H22I1 using the CD-Bonn nucleon-nucleon 
potential [39]. All partial waves with the total angular momen- 
tum of two-nucleon subsystem up to j max -5 were taken into 
account. One can see that the results of Faddeev calculations 
based on CD-Bonn potential even without involving 3NFs are 
in good agreement with all the analyzing powers within the 
achieved experimental accuracy. 

The dashed lines correspond to a relativistic calculation in the 
multiple scattering expansion formalism ll40Tl up to the second- 
order terms of the nucleon-nucleon f-matrix 114 ill with the use 
of the CD-Bonn Q DWF. The parameterization JH of the 
AW f-matr ix ||4rll obtained from the recent phase-shift analysis 
data SP07 [43] allows to avoid the convergence problem due to 
maximal number of partial wave states in the AW system. The 
model i40Tl takes into account the off-energy-shell effects. The 
approach describes reasonably well the angular dependencies 
of Ay over the whole angular range of measurements and A vv at 
backward angles. It fails to reproduce the behavior of A xx . 

The dot-dashed curves correspond to the relativistic calcula- 
tion of the optical potential framework [44] up to the total angu- 
lar momentum / = 39/2. The results are obtained with the use 
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Figure 3: The angular dependence of the dp- elastic scattering cross section 
at Tf -880 MeV. The solid, dashed and dash-dotted curves are the results 
of the Faddeev calculation [22], of the relativistic multiple scattering model 
(40ll46fl and of the optical potential approach |44|| . respectively. The open tri- 
angles and circles correspond to the data obtained at T^*=850 MeV 14711 and 
r'"*=940 MeV (H, respectively. 



Figure 2: The angular dependencies of the vector Ay and tensor analyzing pow- 
ers A y y and A xx foi 
plained in the text. 



ers Ayy and A xx for dp- elastic scattering at rj*=880 MeV. The lines are ex- 



of DWF derived from the dressed bag model of the Moscow- 
Tubingen group 14511 and the on-shell nucleon-nucleon f-matrix 
based on the recent phase-shift analysis data SP07 14311 . These 
calculations reproduce the behavior of A vv only at the angles 
larger than 100° in the cm., while they fail to describe the ana- 
lyzing powers A y and A xx . 



The results of the nonrelativistic Faddeev calculations 12211 . 
the relativistic multiple scattering model [40, 46] and the optical 
potential approach [44] for the dp- elastic scattering differential 
cross section at 7^=880 MeV are shown in Fig. [3]by the solid, 
dashed and dash-dotted curves, respectively. They are com- 
pared with the experimental data obtained at 7^=850 MeV 
1 13] and T' d ab =940 MeV JUt] given by the open triangles and 
circles, respectively. 

The optical potential approach [44] which uses the DWF de- 
rived from the dressed bag model 114511 fails to reproduce the 
cross section data. Both non-relativistic Faddeev [22] and rela- 
tivistic multiple scattering ll46ll calculations describe the exper- 
imental data up to the scattering angles in the cm. of ~70°. The 
Faddeev approach [2211 cannot reproduce the cross section data 
at larger angles, while the relativistic multiple scattering model 
calculations ll46Tl provide better agreement with the data in the 
vicinity of the differential cross section minimum. Note that 
the relativistic effects in the framework of Faddeev calculations 
performed with NN forces are significant only in the region of 



backward angles and they are small in the region of the differ- 
ential cross section minimum 12811 . Therefore, the relativistic 
effects alone cannot explain the difference between the predic- 
tions of these two models. Both approaches cannot reproduce 
the data at backward scattering angles, where the essential con- 
tribution of the A-isobar excitation into the differential cross 
section at high energies was demonstrated l49ll . 

The agreement between the cross section data and theoretical 
calculations at the angles larger than 70° in the cm. could be 
improved, if the three-nucleon forces are taking into account. 
However, theoretical predictions including present 3NF models 
1 2C , |21J or an effective 3NF due to explicit A-isobar excitation 
1 2J underestimate the data at 250 MeV/nucleon OllQil] by up 
to 40%. The inclusion of present 3NFs in the calculations does 
not improve the description of the data on the nucleon analyzing 
power and proton polarization transfer coefficients obtained at 
250 MeV IB [3. 

Therefore, at higher energies, the deviation 
of both spin observables and cross sections from the theoret- 
ical calculations 11221 l27l |28J demonstrates the deficiencies of 
the present 3NF models and relativistic description of the dp 
-elastic scattering process. This might indicate that additional 
short-range 3N forces should be added to the 2jt -exchange type 
forces I28II29I1. 

The relativistic multiple scattering model describes reason- 
ably well the dp- elastic scattering differential cross section and 
vector analyzing power A v up to 140° in the cm. HE- How- 
ever, it reproduces the behavior of the tensor analyzing power 
Ayy only at backward angles, while the angular dependence of 
A xx is not described. The agreement between the experimental 
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data and theoretical description could be improved, if the 3NF 
are included into consideration. 



4. Conclusions 

The vector A y and tensor analyzing powers A vy and A xx for 
dp- elastic scattering have been measured for the first time at 
Internal Target Station at the JINR Nuclotron at 7^=880 MeV 
over the cm. angular range from 60° to 140° corresponding to 
the transverse momenta of ~400-600 MeV/c. 

New results on the tensor analyzing powers indicate strong 
deviations from the predictions of the relativistic phenomeno- 
logical approaches |4^, 44] based on the use of the nucleon- 
nucleon forces only. The non-relativistic Faddeev calculations 
1 22 ] using CD-Bonn 2NF J39ll describe only the angular behav- 
ior of the new data for the analyzing powers. They, however, 
fail to reproduce the differential cross section data obtained in 
earlier experiments 11471 14811 . 

Some deficiencies in the description of the differential cross 
section and the deuteron analyzing powers at 7j* ~880 MeV 
obtained at quite large transverse momenta require the consid- 
eration of additional mechanisms, for instance, 3NFs. Since 



present 3NFs models cannot improve the agreement with the 
data obtained at lower energies lOMl^], new models of 3NFs 
(including short-range part) should be considered. 
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